LANDAU EXPANSION FOR THE KUGEL-KHOMSKII t 2g HAMILTONIAN 
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The Kugel-Khomskii (KK) Hamiltonian for the titanates describes spin and orbital superex- 
change interactions between d ions in an ideal perovskite structure in which the three t2 g orbitals 
are degenerate in energy and electron hopping is constrained by cubic site symmetry. In this pa- 
per we implement a variational approach to mean-field theory in which each site, i, has its own 
n x n single-site density matrix p(i), where n, the number of allowed single-particle states, is 6 (3 
orbital times 2 spin states). The variational free energy from this 35 parameter density matrix is 
shown to exhibit the unusual symmetries noted previously which lead to a wavevector-dependent 
susceptibility for spins in a orbitals which is dispersionless in the g a -direction. Thus, for the cu- 
bic KK model itself, mean-field theory does not provide wavevector 'selection', in agreement with 
rigorous symmetry arguments. We consider the effect of including various perturbations. When 
spin-orbit interactions are introduced, the susceptibility has dispersion in all directions in q-space, 
but the resulting antiferromagnetic mean-field state is degenerate with respect to global rotation of 
the staggered spin, implying that the spin-wave spectrum is gapless. This possibly surprising con- 
clusion is also consistent with rigorous symmetry arguments. When next-nearest-neighbor hopping 
is included, staggered moments of all orbitals appear, but the sum of these moments is zero, yielding 
an exotic state with long-range order without long-range spin order. The effect of a Hund's rule 
coupling of sufficient strength is to produce a state with orbital order. 

PACS numbers: 75.10.-b, 71.27.+a 



I. INTRODUCTION 

High temperature superconductivity 1 and colossal 
magnetoresistance 2 have sparked much recent interest in 
the magnetic properties of transition metal oxides, par- 
ticularly those with orbital degeneracy. 3,4 In many transi- 
tion metal oxides, the d electrons are localized due to the 
very large on-site Coulomb interaction, U . In cubic oxide 
perovskites, the crystal field of the surrounding oxygen 
octahedra splits the d-orbitals into a two-fold degener- 
ate e g and a three-fold degenerate ti g manifold. In most 
cases, these degeneracies are further lifted by a cooper- 
ative Jahn- Teller (JT) distortion, 3 and the low energy 
physics is well described by an effective superexchange 
spin-only model. 5-7 However, some perovskites, such as 
LaTi03, 8 ' 9 do not undergo a significant JT distortion, in 
spite of the orbital degeneracy. 10 In these systems, the ef- 
fective superexchange model must deal with not only the 
spin degrees of freedom but also the degenerate orbital 
degrees of freedom. 3 ' 4,11 The large degeneracy of the re- 
sulting ground states may then yield rich phase diagrams, 
with exotic types of order, involving a strong interplay 
between the spin and orbital sectors. 4,8,9 

In the idealized cubic model for the titanates, there 
is one d electron in the t^g degenerate manifold, which 



contains the wavefunctions \X) = d yz , \Y) = d xz , and 
\Z) = d xy . Following Kugel and Khomskii (KK), 11 one 
starts from a Hubbard model with on-site Coulomb en- 
ergy U and nearest-neighbor (nn) hopping energy t. For 
large U, this model can be reduced to an effective su- 
perexchange model, which involves only nn spin and or- 
bital coupling, with energies of order e = t 2 /U. This 
low energy model has been the basis for several the- 
oretical studies of the titanates. In particular, it has 
been suggested 12 that the KK Hamiltonian gives rise 
to an ordered isotropic spin phase, and that an energy 
gap in the spin excitations can be caused by spin-orbit 
interactions. 13 However, these papers are based on as- 
sumptions and approximations which are hard to assess. 
Recently 14 (this will be referred to as I) we have pre- 
sented rigorous symmetry arguments which show several 
unusual symmetries of the cubic KK Hamiltonian. Per- 
haps the most striking symmetry is the rotational invari- 
ance of the total spin of a orbitals (where a — X,Y, 
or Z) summed over all sites in a plane perpendicular to 
the a-axis. This symmetry implies that in the disordered 
phase the wavevector-dependent spin susceptibility for a 
orbitals, Xa( c l) is dispersionless in the g Q -direction. In 
addition, as discussed in I, this symmetry implies that 
the system does not support long-range spin order at 
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any nonzero temperature. Thus the idealized cubic KK 
model is an inappropriate starting point to describe the 
properties of existing titanate systems. This peculiar ro- 
tational invariance depends on the special symmetry of 
the hopping matrix element and it can be broken by al- 
most any perturbation such as rotation of the oxygen 
octahedra. Here we consider the effect of symmetry- 
breaking perturbations due to a) spin-orbit interactions, 
b) next-nearest-neighbor (nnn) hopping, and c) Hund's 
rule coupling. According to the general symmetry argu- 
ment of I, although long-range order at nonzero tempera- 
ture is possible when spin-orbit interactions are included, 
the system still possesses enough rotation symmetry that 
the excitation spectrum should be gapless. (This con- 
clusion is perhaps surprising because once spin-orbit in- 
teractions are included, the system might be expected 
to distinguish directions relative to those defined by the 
lattice.) This argument would imply that mean-field the- 
ory will produce a state which has a continuous degen- 
eracy associated with global rotation of the spins. The 
purpose of this paper is to implement mean-field theory 
and to interpret the results obtained therefrom in light 
of the general symmetry arguments. We will carry out 
this analysis using the variational properties of the den- 
sity matrix. In a separate paper 15 (which we will refer 
to as III, the present paper being paper II) we will study 
the self-consistent equations of mean-field theory which 
contain information equivalent to what we obtain here, 
but in a form which is better suited to a study of the 
ordered phase. Here our analysis is carried out for the 
cubic KK Hamiltonian with and without the inclusion of 
the symmetry-breaking perturbations mentioned above. 
In the presence of spin-orbit interactions we find that 
the staggered moments of different orbital states are not 
collincar, so that the net spin moment is greatly reduced 
from its spin-only value. The effect of nnn hopping is also 
interesting. Within mean-field theory, this perturbation 
was found to stabilize a state having long-range staggered 
spin order for each orbital state, but the staggered spins 
of the three orbital states add to zero. When only Hund's 
rule coupling is included, mean-field theory predicts sta- 
bilization of long-range spin and orbital order. However, 
elsewhere 16 we show that fluctuations favor spin-only or- 
der. As a result, a state with long-range order of both 
spin and orbital degrees of freedom can only occur when 
the strength of the Hund's rule coupling exceeds some 
critical value which we can not estimate in the present 
formalism. 

Briefly this paper is organized as follows. In Sec. II 
we discuss the KK Hamiltonian and fix the notation we 
will use. In Sec. Ill we discuss the construction of the 
mean-field trial density matrix as the product of single- 
site density matrices, each of which acts on the space of 
six one-electron states of an ion, and whose parametriza- 
tion therefore requires 35 parameters. Here we show that 
the wavevector-dependent spin susceptibilities which di- 
verge as the temperature is lowered through a critical 
value have dispersionless directions, so that unusually 



mean-field theory provides no 'wavevector selection' at 
the mean-field transition. In Sec. IV we discuss the 
Landau expansion at quartic order. In Sec. V we treat 
several lower symmetry perturbations, namely spin-orbit 
interactions, nnn hopping, and Hund's rule coupling. In 
each of these cases 'wavevector selection' leads to the 
usual two-sublattice structure, but the qualitative nature 
of ordering depends on which perturbation is considered. 
In Sec. VI we summarize our work and discuss its impli- 
cations. 



II. THE HAMILTONIAN 

The system we treat is a simple cubic lattice of ions 
with one d electron per ion in a d-band whose five orbital 
states are split into an e g doublet at high energy and a 
tig triplet at low energy. Following the seminal work of 
Kugel and Khomskii 11 (KK), we describe this system by 
a Hubbard Hamiltonian Hh of the form 

Km = e °-c\ a(J c laa + ^ ^ t a0 (i,j)c\ aa c j0a 

iacr (ij) after 

+ U J2J2J2 c\ aa c iarT c\ Prj ,ci 0al , (i) 

i a<f3 erer' 

where c\ aa creates an electron in the orbital labeled a 
in spin state a on site i, e a is the crystal field energy of 
the a orbital, t a p(i,j) is the matrix element for hopping 
between orbital a of site i and orbital (3 of site j, and (ij) 
indicates that the sum is over pairs of nearest neighbor- 
ing sites i and j on a simple cubic lattice. It is convenient 
to refer to the orbital state of an electron as its 'flavor'. 
In this terminology cJ Qcr creates an electron of flavor a 
and ^-component of spin a on site i. Initially we consider 
the case when the Coulomb interaction does not depend 
on which orbitals the electrons are in. In a later sec- 
tion we will consider the effects of Hund's-rule coupling. 
In a cubic crystal field, the crystal-field energy e a splits 
the five orbital d states into a low-energy triplet, whose 
states are d yz = X, d xz = Y, and d xy = Z, and a high 
energy doublet, whose presence is ignored. In this model 
it is assumed that hopping occurs only between nearest 
neighbors and proceeds via superexchange through an in- 
tervening oxygen p orbital, so that the symmetry of the 
hopping matrix is that illustrated in Fig. 1. Thus t a p 
is zero if a ^ j3 and t aa (i,j) = t, except that t aa (i,j) 
vanishes if the bond (ij) is parallel to the a-axis. 11 The 
a-axis is called 17 the inactive axis for hopping between a 
orbitals. When t <C U, KK reduced the above Hubbard 
Hamiltonian to an effective Hamiltonian for the manifold 
of states for which each site has one electron in a ti g or- 
bital state. We will call this low-energy Hamiltonian the 
KK Hamiltonian and it can be regarded as a many-band 
generalization of the Heisenberg Hamiltonian. The KK 
Hamiltonian is often written in terms of spin variables 
to make the analogy with the Heisenberg model more 
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apparent, but for our purposes it is more convenient to 
write the (KK) Hamiltonian in the form 

(ij) P-Y^W) VP 

= e Yl ^Qprn-fpi^QiPiPvti) i ( 2 ) 

(ij) 0-yjt(ij) VP 

where e = t 2 /U and the notation /?7 ^ (ij) indicates that 
in the sum over [3 and 7 neither of these arc allowed to 
be the same as the coordinate direction of the bond (ij) . 
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tive axis for a-hopping). In I we showed that the total 
number of electrons in an a-plane which are in a or- 
bitals is constant. In addition, the total spin vector (as 
well as its z component) summed over all electrons in a 
orbitals in any given a-plane was shown to be a good 
quantum number. The fact that one can rotate the spin 
of all a electrons (these are electrons in a orbitals) in 
any a-plane at no cost in energy implies that there is 
no long-range spin order at any nonzero temperature. 14 
Nevertheless, since experiment 8 shows that LaTi03 does 
exhibit long-range spin order, it must be that spin order- 
ing is caused by some, possibly small, symmetry break- 
ing perturbation, which should be added to the ideal- 
ized KK model. Therefore it is worthwhile investigat- 
ing what form of long-range order results when possi- 
ble symmetry-breaking perturbations are included. Al- 
though the mean-field results we obtain below should 
not be taken quantitatively they may form a qualita- 
tive guide to the type of ordering one might expect for 
more realistic extensions of the above KK model. We 
also noted 6 ' 14 that even when spin-orbit coupling is in- 
cluded, the Hamiltonian has sufficient symmetry that the 
spin-wave spectrum remains gapless. As a result, the gap 
observed 8 in the excitation spectrum of LaTiC>3 can not 
be explained on the basis of the KK Hamiltonian with 
only the spin-orbit interaction as a perturbation. As we 
shall see, these symmetries are realized by the mean- field 
solutions we obtain. 



III. LANDAU EXPANSION AT QUADRATIC 
ORDER 




FIG. 1. A schematic view of the \Z) 
the (indirect) hopping parameter t via intermediate oxygen 
p-orbitals. Positive (negative) regions of wavefunctions are 
represented by dark (light) lobes. In (a) we show that the 
hopping matrix elements between orbitals of different flavors 
are zero. In (b) we show that there is no indirect hopping 
along the z-axis for an electron in the Z-orbital, due to sym- 
metry. 
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pointed out several unusual symme- 



tries of this Hamiltonian. By an a-plane we mean any 
plane perpendicular to the a axis (which is the inac- 



We will develop the Landau expansion of the free en- 
ergy as a multivariable expansion in powers of the full 
set of order parameters necessary to describe the free en- 
ergy arising from the KK Hamiltonian. In this section we 
construct this expansion up to quadratic order in these 
order parameters and thereby analyze the instability of 
the disordered phase relative to arbitrary types of long- 
range order. In later sections we discuss how this picture 
is modified by higher-order terms in the expansion, and 
by the addition of various symmetry-breaking terms into 
the Hamiltonian. 



A. Parametrizing the Density Matrix 

The version of mean-field theory which we will imple- 
ment is based on the variational principle according to 
which the exact free energy is obtained by minimizing 
the free energy functional F(p) as a function of the trial 
density matrix p, which must be Hermitian, have no neg- 
ative eigenvalues, and be normalized by Trp = 1. Here 
the trial free energy is 



F(p) = Tr p(H + kTlnp) 



(3) 
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where the first term is the trial energy and the second is 
— T times the trial entropy, where T is the temperature. 
Mean-field theory is obtained by the ansatz that p is the 
product of single-site density matrices, p(i): 



(4) 



and F(p) is then minimized with respect to the variables 
used to parametrize the density matrix, p(i). Since p(i) 
acts in the space of tig states of one electron, it is a 6 x 6 
dimensional Hcrmitian matrix with unit trace. 

The most general trial density matrix (for site i) can 
be written in the form 



where 



with 



p(i) = li + x(i), 



a/3 pJJ 



Y ap p v (i) = A a f 3 (i)S pv + B af3 (i) ■ a p 



(5) 



(6) 



(7) 



Here a is the Pauli matrix vector, and A a p(i), B^(z), 
B^o(i), and B z a Ji) are 3x3 Hermitian matrices, of which 
the first is traceless. The diagonal terms of the matrix A 
are parametrized for later convenience as 

( _ ai(i) ai(i) . _ Qi(i) _ ai(i) 
xx{ ) ~ V6 V2 ' Vvi ) ~ V6 V2 ' 

A zz (i) = ~A xx {i) - A yy (i), (8) 

such that 

A 2 xx (i) + A 2 yy (i) + A 2 xzz (i) = a\{i) + a 2 2 (i), 
-A 2 xx (i) - A 2 yy (i) + 2A 2 zz (i) = a 2 (i) - a 2 S). (9) 

For any operator 0(i) associated with site i we define 

<0(i)) = Tr[0(i)p] , (10) 

where Tr denotes a trace over the six states \ce,a) of the 
atom at site i with a single ti g electron. Then the di- 
agonal matrix elements of A(i) give the occupations of 
orbital states, 



(11) 



which may be related to the matrix elements of the an- 
gular momentum, L, 

{ Ll{l) f l ) = (N x (i)) = \ + -^ 0l (i) + V2a 2 {i), 

( Lg( 3~ 1 ) = WO) = ^ + ^ai(0 - ^2(i), 
L 2 !'?") - 1 14 

— } = {Nz{l)) = 3 - vT lW ■ (12) 



The off-diagonal matrix elements of A(i) are 

a/3 <J 

= -2«^ A Q(3 (i)e Q(37 , (13) 

a/3 

where e a/ 3 7 is the fully antisymmetric tensor. Similarly, 

{Lp(i)L 7 (i) +L 7 (i)L j8 (i)) = -3^(4^0^ + 4^^} 

= -6^(0 + ^(0]. (14) 

Similarly, the diagonal matrix elements of B 7 (i), B 7 a (i), 
give the thermal expectation value of the 7 component 
of the spin of a-flavor electrons: 

(S ai (i)) = 5>L<^<W = 2B 7 Q «. (15) 

The off-diagonal matrix elements of B 7 (i) are related to 
the order-parameters associated with correlated ordering 
of spins and orbits. 

In general, the density matrix Eq. (5) yields the aver- 
age 



+ ^2( C \ a <j C i0V C \ a >p( A »>[3>{i)5pT + B a >/3>(i) ■ Spr)cif3' T ) 

a'0> 
pr 

= SafjSari/S + Ap^S^ + Bf ja (i) ■ S ^ . (16) 

B. Construction of the Trial Free Energy 

Using the result Eq. (16), we get the trial energy, U, 

(ij) np 

= 2e J2 J2 [A a p{i)Ap a { 3 ) + B a[i {i)-B fia {j)], (17) 

{ij) a[3^{ij) 

where we have used the identity 

y ^ (-B ttl Q 2 • (T pi p 2 )(B a2a3 ■ (Tp 2 p 3 ) 

P2 

^Qia2 ' -B a2 a 3 Sp 1 p 3 ~\~ i&pips ' B aia2 x B a2a3 . (IS) 

Here and below we drop terms independent of the trial 
order-parameters. 

Using Eq. (5) we write the trial entropy as 
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as 
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-TS = kT^2TY(3X 2 (i)-6X 3 (i) + l8X i (i) + ...), (19) with 



where we noted that TrX(i) = 0. The second-order con- 
tribution is found from 

a/3 VP 

a' /3' V'P' 

x 4 a > P > Y a'p'f3>v'( i ) c *f>'v'} = J2J2 Ya Pf^ Y f^ a P^ 

af3 rip 

= J2[A a p(i)Ap a (i) + B afj (i) ■ Bp a (i)}. (20) 

a[3 

At quadratic order the trial free-energy, F — F 2 , is 
thus 

ij a/3 

+B al 3(i)-B 0a {j)}, (21) 

where the inverse susceptibility is given by 

X~l{i,j) = 12fcT<% + 2e 7y (l - 5 ij>a )(l - S ij>p ). (22) 

Here 7^ is unity if sites i and j are nearest neighbors 
and is zero otherwise, and Sij. a is unity if the bond (ij) 
is along the a-direction and is zero otherwise. 

C. Stability Analysis - Wavevector Selection 

We now carry out a stability analysis of the disordered 
phase. At quadratic order in the Landau expansion, pos- 
sible phase transitions from the disordered phase to a 
phase with long-range order are signalled by the diver- 
gence of a susceptibility. Depending on the higher-than- 
quadratic order terms in the Landau expansion, such a 
transition may (or may not) be preempted by a first-order 
(discontinuous) phase transition. So mean-field theory is 
a simple and usually effective way to predict the nature 
of the ordered phase in systems where it may not be 
easy to guess it. To implement the stability analysis we 
diagonalize the inverse susceptibility matrix by going to 
Fourier transformed variables, whose generic definition is 



v i 

F(rO = ^£F(q)e--, 



(23) 



where N is the total number of lattice sites. Then the 
free energy at quadratic order is F 2 = J2q ^2 (q) , where 



F * = ^EE^(q)[^(q)^(-q) 



q a/3 



X"» = 12fcT 
+ 2e£ e-^(l S Rnn , a& )(l - S Rnn ^), 



(25) 



R„ 



where R, m is a vector to a nearest-neighbor site, and a 
is the unit vector in the a-direction. We hence see that 
we have only two kinds of inverse susceptibilities, the one 
for the diagonal elements, namely 



X" Q (q) = 12fcT + 2eJ2 e-^ R ""(l - 5 Rn 



R„ 



12fcT + 2eE4/3 7 (c/3+c 7 ), 
0i 



(26) 



and the second for the off-diagonal matrix elements, 
namely 



= 12fcT + 4eE4/3 7 c 7 , (27) 



where c a = cos(q a a). 

At high temperature all the eigenvalues of the suscep- 
tibility matrix are finite and positive. As the temper- 
ature is reduced, one or more eigenvalues may become 
zero, corresponding to an infinite susceptibility. Usually 
this instability will occur at some value of wavevector 
(or more precisely at the star of some wavevector), and 
this set of wavevectors describes the periodicity of the 
ordered phase near the ordering transition. This phe- 
nomenon is referred to as 'wavevector selection'. In ad- 
dition, and we will later see several examples of this, the 
eigenvector associated with the divergent susceptibility 
contains information on the qualitative nature of the or- 
dering. Here, a central question which the eigenvector ad- 
dresses, is whether the ordering is in the spin sector, the 
orbital sector, or both sectors. If the unstable eigenvector 
is degenerate, one can usually determine the symmetries 
which give rise to Goldstone (gapless) excitations. (We 
will meet this situation in connection with our treatment 
of spin-orbit interactions.) In the present case, we see 
from Eqs. (26) and (27) that the instabilities (where an 
inverse susceptibility vanishes) first appear at kT — 2e/3 
for the diagonal susceptibilities. Consider first the sus- 
ceptibilities for unequal occupations of the three orbital 
states. Making use of Eqs. (8) and (26), we write 



E Xaa(q)- 4 a«(q)- 4 aa(-q) 

a 

[ ai(q) a 2 (q) ] Xn (q) 



ai(-q) 

fl2(-q) 



(28) 



+B a p(q) ■ Bp a (-q)}, 



(24) with the 2x2 susceptibility matrix x n given by 
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X~n(q) = 12fcT + y (hc x + 5c y + 2c,) , 
X^2 2 (q) - 12fcT+2e( Cx + C2/ + 2c z ), 

x^i2(q) = XrT,2i(q) = ^= - c *) • 



B 



(29) 



The instability occurs for both eigenvalues of the in- 
verse susceptibility matrix X~* m (q), but only when the 
wavevector q assumes its antiferromagnetic value Q = 
(it, 7r, it) I a which leads to a two sub-lattice structure (see 
Fig. 2) called the "G" state. The two-fold degeneracy 
is the symmetry associated with rotations in occupation 
number space (N x ), (N v ), and (N z ) with the constraint 
that the sum of these occupation numbers is unity. (At 
quadratic order we do not yet feel the discrete cubic sym- 
metry of the orbital states.) In contrast, the inverse 
spin susceptibility x" 1 OI Eq. (26) has a flat branch 
so that it vanishes for kT = 2e/3 for any value of q a , 
when the two other components of q assume the anti- 
ferromagnetic value 7r I a. This wavevector dependence 
indicates that correlations in the spin susceptibility be- 
come long ranged in an a-plane, but different a-planes 
are completely uncorrelated. Note that beyond the fact 
that there is no wavevector selection in the spin suscepti- 
bility, one has complete rotational invariance in £?2 a (q) 
for the components labeled by 7 independently for each 
orbital labeled a. This result reflects the exact symmetry 
of the Hamiltonian with respect to rotation of the total 
spin in the a orbital summed over all spins in any single 
a-plane. 14 If we restrict attention to the G wavevector 
q = Q, we have complete rotational degeneracy in the 
11 dimensional space consisting of the nine £?^ Q (Q) spin 
order-parameters and the two a„(Q) occupational order- 
parameters. Thus at this level of approximation, we have 
0(11) symmetry! Most of this symmetry only holds at 
quadratic order in mean-field theory. As usual, we expect 
that fourth (and higher) order terms in the Landau ex- 
pansion will generate anisotropies in this 11-dimensional 
space to lower the symmetry to the actual cubic sym- 
metry of the system. As we will see, the anisotropy 
which inhibits the mixing of spin and orbit degrees of 
freedom is not generated by the quartic terms in the free 
energy. Perhaps unexpectedly, as we show elsewhere, 16 
this anisotropy is only generated by fluctuations not ac- 
cessible to mean-field theory. 
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FIG. 2. The two sublattice "G" state which consists of two 
interpenetrating simple cubic lattices on each site of which 
the ions are in a given state, either A or B. 

Dispersionless branches of order-parameter suscepti- 
bilities which lead to an infinite degeneracy of mean- 
field states, have been found in a variety of models, 18-21 
of which perhaps the most celebrated is that in the 
kagome 22 and pyrochlore 23 systems. In almost all cases, 
the dispersionless susceptibility is an artifact of mean- 
field theory and does not represent a true symmetry of 
the full Hamiltonian. In such a case, the continuous de- 
generacy is lifted by fluctuations, which can either be 
thermal fluctuations 24 or quantum fluctuations. 25 Here 
we have a rather unusual case in that the spin suscepti- 
bility has a dispersionless direction (parallel to the inac- 
tive axis) which is the result of an exact true symmetry 
of the quantum Hamiltonian which persists even in the 
presence of thermal and quantum fluctuations. 



IV. LANDAU EXPANSION AT QUARTIC ORDER 

To discuss the nature of the ordered state one may con- 
sider the self-consistent equations for the nonzero order- 
parameters which appear below the ordering tempera- 
ture at kT c = 2e/3 and this is done in III. However, the 
types of possible ordering should also be apparent from 
the form of the anisotropy of the free energy in order- 
parameter space which first occurs in terms in the free 
energy which are quartic in the order-parameters. In 
principle, long-range order is only possible when we add 
to the Hamiltonian terms which destroy the symmetry 
whereby one can rotate arbitrarily planes of spins as- 
sociated with a given orbital flavor. In the next section 
we study several perturbations which stabilize long-range 
order. Although the nature of the ordering depends on 
the perturbation, generically the resulting dispersion due 
to this symmetry-breaking perturbation stabilizes the G 
structure, so that the instabilities arc confined to the 
wavevector q = Q. In this section we implicitly assume 
this scenario. 
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Accordingly, we now evaluate all terms in the free en- 
ergy which involve four powers of the critical variables 
B^ a (Q) and A aa (Q) at the wavevector associated with 
the assumed two sub-lattice, or G, structure. These 
terms arise from two mechanisms. The first contribution, 
which we denote F^ 4 \ arises from "bare" quartic terms 
in Eq. (19). The second type of contribution arises in- 
directly through X 3 (i) in Eq. (19). There we have con- 
tributions to the free energy which involve two critical 
variables and one noncritical variable (evaluated at zero 
wavevector). When the free energy is minimized with 
respect to this noncritical variable, we obtain contribu- 
tions to the free energy which are quartic in the critical 

(3) 

order-parameters and which we denote F 4 . 



A. Bare Quartic Terms, F 4 <4) 

The bare quartic terms are obtained from Eq. (A4), 
by taking into account only diagonal matrix elements of 
the matrices A and B. Since the fourth-order term of 
the entropy is multiplied by 18fcT, [see Eq. (19)], and 
we can safely put here 18fcT=12e, we find that the bare 
quartic terms are given by 



24e EE \< a (i) + Q< a (i)sl(i) + si(i) , (30) 



where we have denoted 

s 2 a (i) = {B x aa (i)f + (J3L«) 2 + (B 7 aa (i)f. 



(31) 



Introducing Fourier transformed variables via Eq. (23) 
we thereby obtain terms quartic in the critical order pa- 
rameters as 



F, 



(4) 



(32) 



where now all order parameters are to be evaluated at 
wavevector Q. Using for the matrix elements of A the 
parametrization Eq. (8), we find 

F[ i] = ^{l2(a? + a 2 2 ) 2 + 48V3aia 2 ( S 2 - s 2 y ) 
+ 48(a 2 + al){sl + 4 + s z ) + 24(4 + 4 + 4) 

-24( a ?- a 2)( s 2 +s 2_ 2s 2)}. (33) 



B. Induced Quartic Terms, F± 

To obtain the terms of this type, we first take from 
Eq. (A2) all the terms having diagonal matrix elements. 
Multiplying them by — 6kT = — 4e [see Eq. (19)], we have 



V 3 = -8e 



+ 3A aa (i)B aa (i) ■ B aa (i) . (34) 



Next we insert here the Fourier transforms. The criti- 
cal variables we treat here are the Fourier components at 
wavevector Q = (tt, 7r, 7r)/a. When the wavevector is Q, 
it will be left implicit. We indicate explicitly only those 
variables taken at zero wavevector. Then V3 is given by 



+ 2B aa (0) ■ B aa A aa , (35) 

where we have used Eq. (31). 

We now eliminate the noncritical variables at zero 
wavevector by minimizing the free energy with respect 
to them. We note that all the noncritical zero wavevec- 
tor variables have the same susceptibility 

X (0) - (12fcT + Se)- 1 = (16c)" 1 , (36) 

and therefore the function to minimize is 

^3 - V 3 + 8e E [A 2 aa (0) + B aa (0) ■ -B aa (0)J . (37) 



The minimization procedure, allowing for the constraint 
J2 a A aa (0) = 0, yields 



-^aa(O) — fjy-B^aAaa, 



A**(o) = ^=(2^ + 24 ■ 
1 



A yy (0) = 



2VN 



2At 



A 2 

A yy 


-4 


- A 2 - 


A 2 

XX 


-4 


~ A zz — 



A zz (0) = -A xx (0) - Ayy(0) 



(38) 



Inserting these values into Eq. (37) yields the contribu- 

(3) 

tion F4 to the free energy 

- J2± V A 2 s 2 - — PTf4 2 +s 2 ) 2 



(•^11 + 4)(A yv + (A yy + Sy)(A zz + S 2 .) 



(39) 



(A zz + S Z )(A XX + s x ) 
which, upon inserting the parametrization (8) becomes 
if) = ^{-12(E 4) 2 + 36 £ 44 ~ 3(a? + a 2 ) 2 

a a<f3 

- 36V3a ia2 (s 2 - s 2 ) - 24 ( a i + a l)(4 + 4 + 4) 
+ 18(a 2 - a 2 )(.s 2 + S 2 -2 S 2 )}. (40) 
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C. Total Fourth-Order Anisotropy 

Adding and F 4 (4) , wc find i<4 as 

F * = ^{ 12 (E s l) 2 - 12 E s «4 + 9(«? + «2) 2 

a a</3 

+ 24(a 2 + a 2 ) ^ s 2 + Uy/Sa^sl - s 2 ] 

a 

-&{a\-al)[s 2 x + sl-2sl]} , (41) 

where all variables are evaluated at wavevector Q. As 
mentioned above, the anisotropy of this form determines 
the nature of the mean-field states of the ideal KK Hamil- 
tonian. We will give a complete analysis of the symmetry 
and consequences of this fourth order anisotropy in paper 
III. Here we will use this form to determine the nature 
of possible ordered states in the presence of symmetry- 
breaking perturbations such as the spin-orbit interaction. 



V. SYMMETRY-BREAKING PERTURBATIONS 

As we have just seen, the idealized KK model con- 
sidered above has sufficient symmetry that there is no 
wavevector selection 26 within mean-field theory and the 
exact symmetry of this model does not support long- 
range order at nonzero temperature. In this section wc 
consider the effects of various additional perturbations 
which are inevitably present, even when there is no dis- 
tortion from perfect cubic symmetry Wc consider in turn 
the effects of a) spin-orbit coupling, b) further neighbor 
hopping, and c) Hund's rule or Coulomb exchange cou- 
pling. Here we do not assume that the long-range order 
only involves the wavevector Q of the G structure. In 
other words our first objective is to see how these vari- 
ous perturbations lead to (if they do) wavevector selec- 
tion and what types of ordering result. 



A. Spin-Orbit Interactions 

We first consider the effect of including spin-orbit in- 
teractions, since these interactions destroy the peculiar 
invariancc with respect to rotating planes of spins of dif- 
ferent orbital flavors independently. Below we see that 
the addition of spin-orbit coupling leads to a wavevector 
selection from the susceptibility, which previously had a 
dispersionless axis in the absence of such a perturbation. 
Indeed, a plausible guess is that the system will select 
the wavevector Q to allow simultaneous condensation of 
spins of the all three orbitals. 

We write the spin-orbit interaction, Vso, as 

^so = A^EE^I^I^a^M/-. ( 42 ) 

i af3~f fJ-v 



where 

(a|L 7 |/3) = -ie aj0 , (43) 

and A is the spin-orbit coupling constant. We now in- 
corporate this perturbation into the mean-field treat- 
ment. The expression for the entropy does not need to 
be changed. The trial energy involves Tr[p(i)Vso] and 
generates a perturbative contribution to the free energy 
which is 

In terms of Fourier transformed variables this is 

SF = 2XN^J2Kp(l = OML»- (45) 

a/37 

Thus the spin-orbit interaction appears as a field act- 
ing on the noncritical order-parameter B aj 3(q = 0), with 

We now calculate the perturbative effect of the spin- 
orbit interaction. Because the perturbation Vso is the 
only term in the Hamiltonian that causes a transition 
from one orbital to another, the leading perturbation to 
the free energy will be of order A 2 . We develop an expan- 
sion at temperatures infinitcsimally below T c = 2e/(3fc) 
in powers of A and {ip}, where {ip} denotes the set of 
variables which, in the absence of spin-orbit coupling, are 
critical at the highest temperature, namely, kT = 2e/3. 
This set includes £?^ Q (q) for q on its "soft line" , which is 
q a arbitrary and the other components equal to ir/a. In 
addition, this set also includes A Qa (Q), namely, ai(Q) 
and a<2, (Q) . The dominant perturbation to the free energy 
will be of order \ 2 ipiipj, where tpi is one of the critical or- 
der parameters. Terms of order X 2 ipi are not allowed, as 
they would cause ordering at all temperatures above T c 
and contributions independent of tpi are of no interest to 
us. So our goal is to calculate all terms of order X 2 tpitpj. 
By modifying the terms quadratic in the critical order 
parameters we will obtain a free energy without a dis- 
persionless branch of the susceptibility, and therefore the 
spin-orbit perturbation will lead to wavevector selection. 

Terms of order X 2 tpitpj in the free energy arise from 
either bare fourth-order terms or indirectly from cubic 
terms which involve one noncritical variable and two crit- 
ical variables. Here we describe these contributions qual- 
itatively. The explicit calculations are given in Appen- 
dices B and C. We first consider contributions arising 
from the third-order terms. Note that the spin-orbit per- 
turbation Vso acts like a "field" in that it couples linearly 
to the order parameter B^(q = 0), as one can sec from 
Eq. (45). Minimization with respect to this order pa- 
rameter yields 

KM = °) = -^ 1/2 (/3|£» = iN^goecfr, (46) 
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where go = A/(6e) and we noted that the non-diagonal 
inverse susceptibility Xafii®) 1S 12e at fcT = 2e/3 [see Eq. 
(27)]. In other words, we have the spatially uniform dis- 
placement, B^g(i) = igoe a f}-y, which is linear in A. Now 
consider third-order terms in the free energy which are 
schematically of the form 



SF 



0)tpiX 



J ■ 



(47) 



where a is a constant, and Xj IS cL noncritical variable, so 
that its susceptibility \j is finite at T c . After minimizing 
with respect to Xj, we obtain a contribution to the free 
energy of order — (l/2)xja 2 (B^ ( g(q = 0)) 2 ip 2 , which is a 
term of order A 2 ^^ (albeit with i = j). This pertur- 
bative contribution to the free energy quadratic in the 



critical variables will be denoted R 



(3) 



Note that these 



cubic terms [see Eq. (47)] are identified as being linear 
in (a) BLj, in (b) a critical order-parameter if>i, such as 

B aa (q a ) (by this we mean B aa evaluated for a wavevec- 
tor on its soft line), or A aa (Q), and in (c) some non- 
critical order-parameter. Terms of order X 2 ipiipj can also 
come from bare fourth order terms which are products 
of two powers of B^g(q = 0) wrt h t wo critical variables 

and these contributions are denoted i 7 ^ 4 ' . All these terms 
will then lead to modifications of the terms in the free 
energy which are quadratic in the critical variables and 
which therefore may lead to wavevector selection within 
the previously dispersionless critical sector. 

We now identify cubic terms in Eq. (19) which are 
of the form written in Eq. (47). There are no nonzero 
cubic terms which are linear in both A and either ai(Q) 
or a2(Q). The allowed cubic terms are analyzed in Ap- 
pendix B and the result for their pcrturbative contribu- 
tion 5P 2 to the free energy from minimizing these cubic 
terms is 



F 2 (3) = -C ^e 2 /37 {^(2 Sa7 (q) Sa7 (-q) 

a/37 la 

+ s aa (q)s aa (-q)j 

+ (s Q7 (Q)WQ) - 2s /3/3 (Q)s Q/3 (Q)) }, (48) 

where Co = 144g 2 e = 4A 2 /e, and we have introduced the 
definition 



s a/3 (q) = B£ Q (q). 



(49) 



In Eq. (48), means that the wavevector is summed 
over the soft line so that = ir/a for \i ^ a and q a 
ranges from — n/a to n/a. In particular the sum over q a 
also includes q = Q. In Appendix C we evaluate the bare 
quartic terms in the free energy which also give a result 
of order X 2 \piipj, and find 

F 2 ] = c o{lJ2 E Wq)W-q) + o?(Q) + 



l E <U (2WQ) WQ) - E «a„(Q)WQ)) }• 



a/37 



(50) 



We now discuss the meaning of these results. One ef- 
fect of the spin-orbit contributions is to couple critical 
spin variables of different orbitals. But this type of cou- 
pling only takes place at the wavevector Q at which spin 
variables for both orbitals are simultaneously critical. So 
we write the sum of all the quadratic perturbations in 
terms of spin variables s Q7 listed above as 

6F2 ^ ^EE[(Et M d Q VWq)W-q)) 
+[M< a) ] F v(Q) s -(Q)]. ( 51 ) 

where Mj"-* is a diagonal matrix and Mi Q ' is an off- 
diagonal matrix. These matrices are 



(a) 



M 



Ml"' = 



4C 
3 

4C 



1 
1 
1 



1 1 

1 -1 
1 -1 



(52) 



where the first row and column refers to s aa and the 
other two refer to sp a , with (i ^ a. The contributions to 

the free energy from M^"' arc independent of wavevec- 
tor and thus do not influence wavevector selection. The 
term in M" selects Q (because the minimum eigenvalue 

of the matrix Mo"' 1 is — 4Co/3, which is negative). In 
addition, the minimum eigenvector determines the linear 
combination of order parameters that is critical. If this 
eigenvector has components (ci,C2,c 2 ), then, for a = x, 
we have 

Sxx(Q) = £, x ci , s yx {Q) = i x c 2 , s zx (Q) = £ x c 2 , (53) 

where £ x is the normal mode amplitude and we adopt 
the normalization cf + 2c?, = 1. Thus, out of the nine 
spin components s a| g(Q) which were simultaneously crit- 
ical in the absence of spin-orbit coupling, we have the 
spin fluctuation corresponding to the three normal-mode 
amplitudes £ x , £ y , and £ z in terms of which we write the 
staggered spin vector for orbital a, s a (Q), as 

Sx(Q) = {Sxx(Q),S X y(Q),S xz (Q)) = (£xCl,€yC2,€ z C2), 
s z(Q) = {€xC2,€yC2,£zCl). (54) 

The total spin at site i is the sum of the spins associated 
with each orbital flavor and is given by the staggered spin 
vector 



S(Q) = +2C2) , 



(55) 
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so that the £'s are proportional to the components of the 
total spin. Now we evaluate the fourth-order free energy 
terms relevant to the spin order-parameters [see Eq. (41)] 
in terms of these critical order parameters 

sf = a ([£ + e y + e z fl4 + $4 + 2c\ c \\ 

where C\ is a constant. In general, a form like this would 
have "cubic" anisotropy in that the vector £ (the total 
spin vector) would preferentially lie along a (1, 1, 1) di- 
rection in order to maximize the negative term in 
However, for the present case, the minimum eigenvector 

of Mo"' 1 is (ci, c 2 , c 2 ) oc (1, — 1, —1). Thus for the present 
case c\ = c 2 ., and the quartic term is isotropic in £ space. 
What this means is that although the spin-orbit interac- 
tion selects the directions for the spin vectors s a of or- 
bital flavor a relative to one another, there is rotational 
invariance when all the s Q 's are rotated together. This 
indicates that relative to the mean-field state there are 
zero frequency excitations which correspond to rotations 
of the staggered spin. Here we find this result at order A 2 . 
More generally, one can establish this rotational invari- 
ance to all orders in A and without assuming the validity 
of mean-field theory. 14 ' 6 

Note that the spin state induced by spin-orbit coupling 
(with ci = — c 2 ) does not have the spins of the individ- 
ual orbitals, s a , parallel to one another and thus the net 
spin, S, is greatly reduced by this effect. Explicitly, when 
ci = — C2, we have 

S 2 = + $ + 

= s 2 (Q) = s 2 y (Q) = s 2 (Q) = (C + e y + aA (57) 

This means that the total spin squared is 1/3 of what it 
would be if the s a were parallel to one another. 

It remains to check that the variables afc(Q) are less 
critical than s a7 (Q). The results given in Eq. (C2) of 
Appendix C show a positive shift in the free energy as- 
sociated with the variables afc(Q), whereas the spin vari- 
ables have a negative shift in free energy due to spin-orbit 
interactions. We therefore conclude that in the presence 
of spin-orbit interactions, mean-field theory does give 
wavevector selection and one has the usual two-sublatticc 
antiferromagnet, but with a greatly reduced spin magni- 
tude. It is interesting to note that 8 LaTiC>3 has a zero 
point moment which is about 45% of the value of the 
spin were fully aligned. This zero-point spin reduction is 
much larger than would be expected for a conventional 
spin 1/2 Heisenberg system in three spatial dimensions. 
It is possible that spin-orbit interactions might partially 
explain this anomalous spin reduction. 

B. Further Neighbor Hopping 

We now consider the effect of adding nnn hopping to 
the Hubbard model of Eq. (1). For a perfectly cu- 



bic system, this hopping process comes from the next- 
to-shortest exchange path between magnetic ions, as is 
shown in Fig. 3. We write the perturbation V to the 
Hubbard Hamiltonian due to these processes as 

V = t'J2^(i,j)V ij: (58) 

a 

where t' is the effective hopping matrix clement connect- 
ing next-nearest neighbors, a is summed over coordinate 
directions x, y, and z, J a {i,j) is unity if sites i and j are 
next-nearest neighbors in the same a-plane and is zero 
otherwise, and 

Vi i = e *f35 c h* c js* ■ (59) 

a j36 

Here a is in the direction normal to the plane containing 
spins i and j, and e^g S restricts the sum over (3 and 5 to 
the two ways of assigning indices so that a, (3, and 5 are 
all different. Note that the paths from if} to j5 and from 
iS to j[3 use alternate paths of the square plaquette con- 
necting i and j. Notice that the processes which couple 
nearest neighbors cancel by symmetry (see Fig. 4), so 
that the effect of hopping between magnetic ions via two 
intervening oxygen ions involves only nnn hopping. This 
generates a perturbation to the KK Hamiltonian (which 
describes the low-energy manifold) of the form 

Vkk = -e'^^7a(i.i) 

a ij 

06a 136a 

where e' = (t') 2 /U and U is the on-site Coulomb energy. 
This may be written as 

a ij 

where, apart from a term which is a constant in the low- 
energy manifold, we have for a = x 

V x (i,j) — ^""^ (clva c 'i-zy c ivn c i za ' + c \ya c iyn c \z"n c i za 

arj 

Jr 4za C izV c \yn C iy cr 4.za C im4zri C 3V<^ ' (^2) 

and similarly for y and z. 
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FIG. 3. Hopping between different orbitals on next nearest 
neighboring (nnn) Ti ions when hopping between neighboring 
oxygen p orbitals is allowed. The hopping matrix element is 
the product of matrix elements to hop from a Ti ion in a d yz 
state to an ion in a p y state, then to an adjacent O ion also 
in a p y state, and finally to a nnn Ti ion in a d xy state. 




FIG. 4. Hopping between different orbitals on near- 
est-neighboring Ti ions when hopping between neighboring 
oxygen p orbitals is allowed. The matrix elements for the 
two channels to hop from d yz to d xy have opposite signs, so 
that the total matrix element (summed over the two channels) 
is zero, as one would deduce from symmetry considerations. 
Thus the only processes involving two nearest neighboring 
oxygen ions are processes like those shown in Fig. 3 between 
nnn Ti ions. 



The details of the mean-field treatment of this pertur- 
bation is given in Appendix D. Here we summarize the 
major analytic results obtained there for the wavevector- 
dependent spin susceptibility at the critical wavevector, 
Q, Xckx;/3<x' (Q) = Xap{Q)S<T,<r', where a and (3 are or- 
bital indices and a and a' are spin indices. The result of 
Appendix D is that 



12fcT-8e 8e' 
8e' YlkT - 8e 
8e' 8e' 



8e' 
12fcT-8e 



The minimum eigenvalue is 



This gives 



A = 12fcT-8e-8e'. 



kT c = 2(e + OA 



(63) 



(64) 



(65) 



By considering the eigenvectors and the effect of the 
fourth order terms, the analysis of Appendix D shows 
that nnn hopping does stabilize a Q antifcrromagnetic 
structure, but the resulting 120° state has zero net stag- 
gered spin. In addition, as before, there is a degeneracy 
between the spin-only states we have just described, and 
a state involving orbital order. As shown in III, fluctu- 
ations remove this degeneracy, so that we may consider 
only the mean-field solutions for spin-only states. Such a 
magnetic structure for which the local moment (summed 
over all flavors) vanishes, will be rather difficult to detect 
experimentally. 

It is instructive to argue for the above results with- 
out actually performing the detailed calculations of Ap- 
pendix D. We expect the effect of indirect exchange be- 
tween nnn's to induce an antiferromagnetic interaction 
between the spins of different orbital flavors of nnn's. 
Note that the wavevector Q describes a two sub-lattice 
structure in which nnn's are on the same sub-lattice. Ac- 
cordingly, as far as mean-field theory is concerned, an 
nnn interaction between different flavors is equivalent to 
an antiferromagnetic interaction between spins of differ- 
ent flavors on the same site. So the spins of the three 
orbital flavors form the same structure as a triangular 
lattice antifcrromagnet, namely the spins of the three 
different orbital flavors are equal in magnitude and all lie 
in a single plane with orientations 120° apart. This state 
still has global rotational invariance, but also, as does 
the triangular lattice antifcrromagnet, it has degeneracy 
with respect to rotation of the spins of two flavors about 
the axis of the spin of the third flavor. 



C. Hund's Rule Coupling 

We now consider the effect of Hund's rule coupling. 
Our aim is to see how this perturbation selects an or- 
dered phase from among those phases which would first 
become critical in the absence of this perturbation as the 
temperature is reduced. To leading order in rj = Jh/U, 
where Jh is the Hund's rule coupling constant (which 
is positive in real systems), as discussed in Appendix E, 
this perturbation reads 28 



sn 



KK 



*>£ E £(< 

(ij) 017 t (ij) aa ' 



,t 
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C i- 1 cr' C iflcr C j' 1 cr C jl3a' + C ! 7 <r C i/3o- c ] fja' C jjcr' 

), (66) 

where e = t 2 /U, as before. 29 To see the effect of this 
perturbation within mean-field theory, we calculate its 
average (see Appendix E for details). Confining to aver- 
ages which are critical when 77 = 0, (i.e., A aa and B aa ), 
the result of Appendix E is 

(SH K K)=ev^ £ (lOApp{i)A P p{j) 

- lOApp{i)A^{ 3 ) + 2B pp {i) ■ B [ip (j) 
-2Bp (i).B^{j)). (67) 

Using Eqs. (8) and (49) to write the order parameters 
in terms of the a^'s and the s Q7 's, this contributes a per- 
turbation to the free energy given by 

SF ^E^X-^^qM-q) 

k,l 

+ \Y, 5 [xjH^llp *a 7 (q)*|8 7 (-q), (68) 



a/37 



where 

%^(q) 

= -20er? 



3 {2C X -\- 2Cy C^) ^7^{Cx C y) 



^( c x c y) 



-c 2 



and 



c z 
c z 



Ok 





(69) 



(70) 



If the minimum eigenvalue of Sx^ 1 at wavevector Q is 
negative, then the instability temperature for the as- 
sociated order parameter is raised by the perturbation 
and vice versa. Note that at wavevector Q, c x — c y — 
c z = — 1 the eigenvalues of 5 [x7 x (q)] are 877c, — 47/e, and 
—A-qe. On the other hand, the eigenvalues of 6 [x^ 1 (q)] 
are both — 20r/e. From this result we conclude that 
Hund's rule coupling favors antiferromagnetic orbital or- 
dering, as described by the order parameters ai(Q) and 
fl 2(Q)- Since the mean- field temperature for spin and 
orbital ordering were degenerate for 77 = 0, we conclude 
that within mean-field theory the addition of an infinites- 
imal Hund's rule coupling gives rise to an ordering tran- 
sition in which the ordered state shows long-range anti- 
ferromagnetic orbital order, characterized by the order- 
parameters ai(Q) and a 2 (Q). However, since we have 
shown elsewhere 16 that for the bare KK model, fluctu- 
ations stabilize the spin-only states relative to orbital 



states, we conclude that when fluctuations are taken into 
account, it will take a finite amount of Hund's rule cou- 
pling to bring about orbital ordering. For spin ordering 
the mean-field state is degenerate with respect to an arbi- 
trary rotation. This is reflected by the fact that the term 
which is fourth order in the spin components is isotropic. 

We now discuss the anisotropy in the mean-field solu- 
tion for orbital order. We want to determine the form the 
free energy assumes in terms of the Fourier-transformed 
variables ai(Q) anda2(Q). Wavevector conservation dic- 
tates that we can have only products involving an even 
number of these variables. If we write ai(Q) = acos#Q 
and 02 (Q) = asin#Q, then we show in Appendix F 
that the contribution to the free energy of order a 4 is 
independent of 8q, but the term of order a 6 is of the 
form SF = a 6 [C + C 6 cos(66» Q + 0)]. This form indi- 
cates an anisotropy, so that the mean-field solution is 
not subject to a rotational degeneracy in a\-a2 space. 
If Cq is positive and = 0, these minima come from 
the six angles that are equivalent to 0q = n/2 + mr/3. 
For 9q = it/2, a\ — and we have ordering involving 
only a 2 , so that (N z ) = 1/3, (N x ) = 1/3 + V2a 2 (i) and 
(N y ) = 1/3 — \/2«2 ) ■ The six minima of cos(6#q) corre- 
spond to the six permutations of coordinate labels which 
give equivalent ordering under cubic symmetry. Some- 
what different states occur for Ce negative, but different 
solutions reproduce the cubic symmetry operations. 



D. Spin-Orbit Interactions and Hund's Rule 
Coupling 

Here we briefly consider the case when we include the 
effects of both spin-orbit and Hund's rule coupling. We 
consider the instabilities at wavevector Q. In this case 
we construct the spin susceptibility (Q) [defined as in 
Eq. (68)]. For the present case we may use our previous 
calculations in Eqs. (52) and (68) to write 



X7 1 (7) 



\ + x y y 
y X + x z 
y z X Q + x 



(71) 



where the first row and column refer to s 77 and the other 
two rows and columns refer to s 01 with (3 7^ 7 and 



4 4 
x = -77C0, y = -C + Aer), z 



--C + 4e V . (72) 



Similarly the orbital susceptibility (also at wavevector Q) 
is given by 



\q + w 
A + w 



where 



w = 2C - 20efy . 



(73) 



(74) 
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In the above Co = 4A 2 /e must be positive, Ao = \2kT + 
8e, and 77 = Jh/U is normally positive, although we may 
draw a phase diagram incorporating the possibility that 
77 is negative. 

As we have seen, with only spin-orbit interactions we 
get a spin state which has a rotational degeneracy, and 
with only Hund's rule interactions, the ordered phase 
has orbital rather than spin ordering. When both in- 
teractions are present, there is a competition between 
these two types of ordering. To study this competition 
we need to compare the minimum eigenvalue of the two 
susceptibility matrices given above. For the inverse spin 
susceptibility matrix y > z, in which case the minimum 
eigenvalue is 

A- = A + x + (z/2) - vW+V • (75) 

On dimensional grounds, we expect that for Co < Trje, 
where t is a constant, Hund's rule coupling will domi- 
nate and will lead to orbital ordering. Indeed after some 
algebra we find this condition with r w 2.7. This may be 
written as 77 > and A < T'e^/rj, where r' — y/r/2 w 0.82 



n 




FIG. 5. The mean-field phase diagram as a function of the 
spin-orbit coupling constant A and the Hund's rule coupling 
constant r\ = Jh/U (which is normally positive). In the 
"spin-only" phase for 77 =fc 0, the staggered moment orients 
along a (1, 1, 1) direction, but the staggered spin moments of 
different orbital states are not collinear, thus reducing the net 
staggered spin. For 77 = 0, the mean-field state has rotational 
degeneracy, so no easy direction of staggered magnetization 
is selected and the excitation spectrum is gap less. In the or- 
bital phase one has the six-fold anisotropy associated with the 
equivalent choices for differently populating orbital levels in 
cubic symmetry, as is discussed in the text. 

which gives rise to the phase diagram shown in Fig. 5. 
This phase diagram is not quite the same as that found 
in Ref. 28 for zero temperature. When we have spin or- 
dering, we may analyze the fourth-order terms, as is done 
in Eq. (56). That analysis shows that unless the min- 
imum eigenvector has components of equal magnitude, 
the anisotropy favors spin ordering along a (1,1,1) direc- 



tion. The condition that the eigenvector be (— 1, 1, 1) is 
that y + z = 0. This can only happen when 77 = 0. Then 
we have isotropy and the mean-field state exhibits rota- 
tional degeneracy. Otherwise, when 77 ^ 0, the fourth- 
order terms give rise to an anisotropy that orients the 
staggered spin along a (1, 1, 1) direction. We should also 
remind the reader that fluctuations favor the spin-only 
state, so that the phase boundary shown in Fig. 5 will be 
shifted by fluctuations to larger positive 77. In the regime 
of orbital ordering, we indicate in Appendix F the exis- 
tence of a six-fold anisotropy in the variables ai(Q) and 
a 2(Q), such that the six equivalent minima correspond 
to the six possible states which are obtained by choosing 
N a = 1/3 for one coordinate a, and then occupying the 
two other orbitals with probability 1/3 ± A. 



VI. DISCUSSION AND SUMMARY 

The cubic KK model has some very unusual and inter- 
esting symmetries which cause mean-field theory to have 
some unusual features. In particular, for the simplest KK 
Hamiltonian, we found that mean-field theory leads to 
criticality for the wavevector-dependent spin susceptibil- 
ity associated with orbital a which is dispersionless along 
the q a direction of wavevector. This result is consistent 
with the previous observation 14 that the Hamiltonian is 
invariant against an arbitrary rotation of the total spin 
in the orbital a summed over all spins in any single plane 
perpendicular to the a axis. This 'soft mode' behavior 
prevents the development of long-range spin order at any 
nonzero temperature, 14 even though the system is a three 
dimensional one. 

Any perturbation which destroys this peculiar sym- 
metry will enable the system to develop long-range spin 
order. In particular, we investigate the role of a) spin- 
orbit interactions, b) second-neighbor hopping, and c) 
Hund's rule coupling in stabilizing long-range spin order. 
In the presence of spin-orbit interaction we find wavevec- 
tor selection (because now the spin of different orbitals 
can not be freely rotated relative to one another) into 
a two-sublattice antiferromagnetic state with a greatly 
reduced spin magnitude. Since experiment shows such 
a reduction, 8 this mechanism may be operative to some 
extent. However, as noted previously, 14 the excitation 
spectrum does not have a gap until further perturbations 
are also included. The mean- field solution is consistent 
with this conclusion, because the mean-field state which 
minimizes the trial free energy is degenerate with respect 
to a global rotation of the staggered spin. 

The ordered state which results when nnn hopping is 
added to the bare KK Hamiltonian is quite unusual. In 
this state, each orbital flavor has a staggered spin mo- 
ment, but these three staggered spin moments form a 
120° degree state such that the total staggered spin mo- 
ment (summed over the three orbital states) is zero! It 
is not immediately obvious how such long-range order 
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would be observed. Finally, we show that when the bare 
KK Hamiltonian is perturbed by the addition of only 
Hund's rule coupling, the resulting ordered state may 
exhibit long-range antiferromagnetic orbital order. 

One caveat concerning our result should be mentioned. 
All our results are based on a stability analysis of the 
disordered phase. If the ordering transition is a discon- 
tinuous one, our results might not reveal such a transi- 
tion. In III we will present results for the temperature- 
dependence of the various mean-field solutions. Further 
analysis of the ordered phase is needed to obtain a phase 
diagram at T = 0, as is done in Ref. 28. 

It should be emphasized again that all the results in 
this paper are based on the assumption that nearest- 
neighbor bonds along an axis a are 'inactive', namely 
that there is no direct hopping between a orbitals along 
such bonds. Even within cubic symmetry, such hopping 
could still exist, alas with a very small hopping energy 
t". However, as soon as we add such terms, the ver- 
tical bond in Fig. lb becomes active, and Eqs. (26) and 
(27) have the additional contributions &Xaa = 2e"c Q and 
Ax~p = 2e"(c Q + eg), with e" = t" 2 /U. This introduces 

dispersion in all directions, and select order at q = Q. 
Distortions away from the cubic structure can enhace t" , 
and stabilize such order even further. 



One general conclusion from our work is that it is not 
safe to associate properties of real experimental systems 
with properties of a model Hamiltonian unless one is ab- 
solutely sure that the real system is a realization (at least 
in all important aspects) of the model Hamiltonian. Here 
the ideal cubic KK Hamiltonian has properties which are 
quite different from those observed for systems it suppos- 
edly describes. What this means is that it will be nec- 
essary to take into account effects that one might have 
been tempted to ignore in order to identify a model that is 
truly appropriate for experimentally realizable systems. 
Alternatively, perhaps our work will inspire experimen- 
talists to find systems that are as close as possible to that 
of the ideal cubic KK Hamiltonian treated here. Such 
systems would have quite striking and anomalous prop- 
erties. 
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APPENDIX A: HIGHER-ORDER TERMS IN THE FREE-ENERGY 



Here we employ Eqs. (5), (6), and (7) in conjunction with Eq. (19), to derive general expressions for the cubic and 
quartic terms of the free energy. 

The 'bare' cubic terms in the free-energy arise from Tr[X 3 ]. We find 



C i03V3 



Tr[X (i)] — ^""^ Tr \c\ a] pj Y aiPl p iril {i)cji3 lVl cl ao ^ Y a2P2/ 3 2ri2 (i)ci0 2rl2 c\ a3p3 Y a3 p 3 p 3ri3 (i)cj i 

|^aia 2 W^Pll»2 + B ai a 2 (i) ■ (T pi p 2 A a2 a 3 (i)Sp 2 p 3 + Ba 2 a 3 (i) ■ (T p 2 p 3 A a3ai (i)6p 3 p 1 + B aaai (i) ■ Cp 3 px ■ (Al) 



Making use of the identity Eq. (18), this becomes 

Tr[A (i)] = 2 ^ ^ -A aia2 (i)A a2a3 (i)A a3ai (i) -j- 3-A aia2 (i)B a2 g 3 (i) •JB Q3Ql (i) 

+i(B aiCC2 (i) x B a2a3 {i)) ■ B a3ai (i)|. 
The 'bare' quartic terms in the free-energy arise from Tr[A 4 ]. We find 

Tr[A 4 00] = ^2 ^^[4 ai p 1 Y c*iPiPim( i ) c iPim4 a2 p2 Y c'2P2P2V2( i ) c i02 



J2V12 



CtiPi PiV: 

XC ia 3 p 3 Y a 3 p 3 3 r] 3 (i) c il3 3 V3 C ia4P4 Y a4P4l3 4 ri4 (i)Cif3 4rl4 



— ^ ^ ^ ] j^aia2 M^plP2 ~l~ B aia2 (l) ■ <J p lP2 A a20l3 (i)Sp 2 p 3 + B a2a3 (l) ■ (Tp 2t , 



Oti p 
X 



j^4a 3 a4 {i)5p 3 p4 + -B Q3Q4 (z) ± Vp 3 p i A a2iai {i)8p i p 1 + B aiai [i) • <jp 



(A2) 



(A3) 
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Again using the identity Eq. (18), this becomes 

Tr[X 4 (i)} =2^| (A aia2 (i)A a2a , 3 (i) + B aiCC2 (i) ■ B a20l3 (if) (A a30l4 {i)A a4ai (i) + B a3Qli (i) ■ B aiai (if) 

~\~ (^A aiCt2 (i)B a20l3 (i) + B aiOC2 (i)A a20l3 (i) + iB aiCt2 (2) x B a2Ct3 (if) 

■ ( A Q3Ct4 (i)B aiai (i) + B a3Ct4 (i)A ai(Xl (i) + iB a3ai (i) x B ai0ll («))}• (A4) 



APPENDIX B: CUBIC FREE-ENERGY TERMS 

Referring to Eq. (A2), the relevant terms for our purpose come from the second and the third terms there. Working 
in Fourier space we hence have 

SF = 7= X! [ 3A «i«2(qi)^a 2Q3 (q2) • Ba 3 ai(-qi - Q2) 



N 

qiq 2 aia2<J3 

(qi) x B a2 „ 3 (q 2 ) • B a3Cll (-qi - q 2 ) . (Bl) 
When one of the quantities B here acts as the spatially uniform field [see Eq. (46)] , this expression becomes 
SF = -= ^ ' ^ ' 3A QlQ2 (q) B Q2Q3 ■ B a;iai (— q) + 3A QlQ2 (q)_B Q2Q3 (— q) • B a3ai 

v q aia2«3 

+ iB aia2 x B a20t3 (q) • B a3Ctl ( q) 4- iB aia2 (q) x B a2a3 ■ B a3Ctl ( q) 

(q) x -B Q2Q , 3 (-q) • JB Q3Ql , (B2) 

where B which does not depend on q is the uniform field. 

We first consider the terms involving the A's. The relevant contributions come from 013 = ct\ [the first term in Eq. 
(B2))] and a 3 = a 2 [the second term there]. Hence we find 



SFa = ^ 2'A»/j(q)S/Ja • (B aa (-q) + %(-q)) , (B3) 



where X)a/3 denotes that a^)3. When we minimize F 2 + SFa with respect to A a p(q), and use Eqs. (27) and (46), 
we get the contribution 

5 F A = -72gleJ2J2 e lpi\- s <*M) + Wq)HW-q) + W-q)P + cos(g 7 a)]-\ (B4) 

q af3j 

where we have defined 

SQ7 (q) = BIM)- (B5) 

Also, since we are interested in the free energy to quadratic order in the order parameters, we have set kT = 2e/3. 

In this result we want to keep only contributions which involve the critical variables. For s Q7 (q) this means that we 
sum over q's such that qp = n/a, for /3 7^ a. Thus for each s aj the wavevector sum is a sum over the component q a , 
with the other components of q equal to n/a. We denote this type of sum by ^2 qa - Furthermore for a term involving 
components s aj and with different orbitals a and (3, this sum reduces to the single wavevector Q = (Tr,n,n)/a. 
So 

sf a = -i^Y,^{Y: tv^t + »«7(Q)WQ)}- (B6) 

a/37 9c ' 

Here we will set [2 + cos(g 7 a)] = 1 because for s ai (with a ^ 7) wc must have q 1 — n/a. This term favors ordering at 
wavevector Q with s a (Q) collinear with s / a(Q), where s Q (Q) is a vector with components [s ax (Q) , s ay (Q) , s az (Q)]. 
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Next we consider the contribution coming from the terms with three B's in Eq. (B2). Here we put one of the 
q-dependent Bs to be diagonal in the orbital indices, to obtain 



SFb = - l ^M £ £' E W^CqX^-q) - B^(-q))B- . (B7) 



,:-!(i/,T. 

q a/3 Qi/3i7 

Eliminating the noncritical B^(q) variables by minimizing F 2 + SFb with respect to them, we get 



SF B = -l296(g kT) 2 J2' ]T X^(q)[s w (q) - s^(q)][(s^(-q) - « a/ j(-q)], (B8) 

a/3 q 

where x is given in Eq. (27), and we have used the definition (B5). As before we set kT — 2e/3 and separate the 
sums to be only over critical wavevectors for each orbital spin vector, in which case we have 

SF B = -144^e^e^ 7 {^[s aa (q)s aQ (-q) + s a/ 3(q)s a/3 (-q)]-2s aQ (Q)s /3a (Q)} . (B9) 

a/37 la 

Here we noted that Xa/3(q) = Xa/3(Q) = l/(4e) because this component of \ depends on q 7 which is always Tr/a in 
the summation over wavevector. 

In summary the total contribution to the quadratic free energy at order A 2 is 

K, (3) - SF A + SF B = -Co £ 4^ (*a 7 (Q)WQ) - 2*ftj(Q)*a/j(Q)) 

a/37 

+ ( S «7(q) S "7 (~q) + S »« (q) s aa (~q) + Sa/3(q)Sa/3 (~q)]) } , (BIO) 

where we set kT = 2e/3 and C = 144.g 2 e. 

APPENDIX C: QUARTIC TERMS IN THE FREE ENERGY 

Now we look at fourth order terms. These involve two critical order parameters and two powers of A. Therefore, 
we pick from Eq. (A4) all terms involving at least two powers of B. Since two of the factors B in each term have to 
be B a p = —Bp a , with a^ft [see Eq. (46)], we see that the terms involving a single power of A vanish. Thus we 
have to consider the expression 

36fer ^ ^ ^ ^ ^4A QlQ2 A a2(X3 B a3CX4 • B a ^ ai -\- 2A aia2 A a3a4 B a2CX3 • B a ^ ai 

a 3 a 4 

~t~ (B aia . 2 • B a20t3 ) (i? a3Q4 • B aAai ) (B aiQ , 2 x B a20l3 ^j • (B a3a4 x B a ^ ai )^ , (CI) 

where A and S are functions of the site index i. The first two members of Eq. (CI) are calculated for the case in 
which the A's are critical, and the B s are given by Eq. (46). Denoting their contribution to the self-energy by 5F^\ 
we find 



<5F 2 (1) = 36fcT^^(4A 2 tQ («) + 2A aa {i)A^(i))B a0 ■ B Pa 

i a/3 

= 36fcT 5o 2 ]Te^ 7 ((4AL« + 2A aa {i)A P0 {i)) = 216fcT 9 2 ]>> 2 (») + a 2 «), (C2) 

a/37 * 

where in the last step we have used Eq. (8). 

(2) 

The contribution of the remaining two members of Eq. (CI) is denoted SF^ . Here we have to take two of the JB's 
as critical, while the other two are given by Eq. (46). To shorten notations, we denote here the critical B as B(i), 
while the non-critical one is simply written as B. We have 
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5F 2 (2) = 72fcT^^[(B aa (i) ■ B aa (i)){B afj ■ B fja ) + {B aa (i) ■ B a[j )(B aa {i) ■ B 0a ) 



i a (3 



+ (B aa (i) ■ B a p){B f 3f 3 {i) ■ Bpa) - (B aa (i) x B aP ) ■ {B 0p {i) x Bp a ) + (B aa (i) x B af3 ) ■ (B aa (i) x B 0a )\ . (C3) 
Making again use of Eq. (46), this expression becomes 



(C4) 



Transforming to Fourier space, noting that only the first term here contains q while in the other two we must 
necessarily have q = Q, (because they involve simultaneous criticality of two flavors), we obtain 

5F^ 2) = 72kTg 2 Q [2 ]T £ s Q/3 (q) Sa/3 (-q) + ]T ^^(Q^Q) - ]T s a „(Q)WQ))] > ( C5 ) 

q a/3 Q/37 v 

where we have used the definition Eq. (B5). The total contribution to the free energy from quartic terms is then 

fP=8fP+8fV>. (C6) 



APPENDIX D: MEAN-FIELD THEORY FOR NNN HOPPING 

Starting from Eq. (60), we may write the perturbation due to next-nearest-neighbors in the form 

VkK = e ' E E j) e »PS [4pp C ^v c ]^ S P + C Up C i^c]s V C jSp 

ij afSS PV 



(Dl) 



Within our mean-field theory, the averages are taken separately on the operators belonging to the site i, and those 
belonging to site j. The required averages are then given in Eq. (16). The following contribution to the trial energy 
U is then 



{V KK ) = 2e'Y J Y.^j) e 

ij af3& 



2 

a/38 



A 50 (i)A S p{j) + B sp {i) ■ B S0 (j) + App(i)A ss (j) + B Pf} (i) ■ B ss {j) . (D2) 



Transforming to Fourier space, noting that each site has four next-nearest neighbors in each a-plane, we obtain 
(Vkk) = 8e' E e ipscpc5 [^(q)^(-q) + %(q) ■ Bf/j(-q) + A w (q)A M (-q) + %(q) • B M (-q)] , (D3) 

q a/3<5 

where eg = 003(5,35). The result Eq. (D3) is now added to Eq. (24), in order to obtain the modifications in the 
susceptibility tensor. Specifying to the diagonal order-parameters A aa and B aa , the susceptibility tensor becomes 
[seeEq. (26)] 



X(q)" 1 



12kT + 4e(c y + c z ) 8e'c x c y 



8e'c x c z 



8e'c x c y YlkT + Ae(c x + c z ) 8e'c y c z 

8e'c x c z 8e'c y c z YlkT + 4e(c x + c y ) 



(D4) 



Now we look at the most critical wavevector, which here is Q. There we have 



x(Q)- 1 



12kT-8e 8e' 
8e' 12kT-8e 



8e' 



8e' 



8e' 12fcT-8e 



(D5) 



We begin with the analysis of the susceptibility tensor of the spin order parameters, which are given by the elements 
of B aa . Then we can use the matrix (D5). The minimum eigenvalue is 



A= 12fcT-8e-8e', 



(D6) 
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which gives 



kT c = 2e/3 + 2e'/3. (D7) 
Correspondingly, there are two degenerate eigenvectors: 

|1> = (0,1,-1)^/2, |2) = (2,-1,-1)^/6. (D8) 



To avoid confusion between orbital and spin labels, we will here denote the orbital states x, y, and z by a, b, and c. 
Then in terms of normal mode vector £ and p we have the orbital spin vectors as 

Sa = -7T' Sb = 7E p + T2^ Sc = 7r~7^' (D9) 

with 

<-3« , + ^ + >-. + (D10 » 

Evaluating the fourth-order free energy [see Eq. (41)] relevant to the spin-order parameters, we find 

What we see is that the fourth-order term does not select a particular direction for order. We have three angles which 
describe the degenerate manifold. For a given value of £ 2 + p 2 , we optimize the term (£ x p) 2 by taking |£| = \p\ and 
making £ perpendicular to p. So, it takes two angles to specify £ (given that its length is fixed) and then we have 
one angle to specify p, given that \p\ — |£| and it is perpendicular to £. We now discuss what this choice of order 
parameters means for the spin vectors. First note that 

s 2 = s 2 = s 2 = 2£ 2 /3 . (D12) 

Also we see that the three orbital spin vectors obey 

s a • s 6 = s a • s c = s 6 • s c = -£ 2 /3 . (D13) 

The three vectors each make a 120° angle with each other and must therefore lie in a single plane. We can fix, say, 
s a . This accounts for two angles. Then the other two spin vectors require another angle to tell which plane they lie 
in. Note that there is zero net staggered moment. There is long-range spin order, but not of any simple type. 

Next we analyze the susceptibility tensor of the occupation order parameters, which are given by the elements of 
A aa . Since the matrix A aa is traceless, we use the parametrization Eq. (8) to obtain from Eq. (D5) the 2x2 matrix 



x^(q) 



12kT + -^(5c x + 5c y + 2c z ) + ^-{c x c y — 2c y c z 2c z c x ) "^j( c f — Cx ^ ~^% Cz ^ c y ~ Cx ^ 

^§( c j/ _ c x) + ^c z (c y - c x ) \2kT + 2e(c x + c y + 2c z ) - &e'c x Cy 



(D14) 



This gives a minimum eigenvalue identical to that of Eq. (D6), which yields the same instability temperature as 
for the spin-only states. However, in the absence of second-neighbor coupling, the spin-only states are favored by 
fluctuations, 16 so that choice should be maintained for infinitesimal next-nearest neighbor hopping. (The situation 
could change when the next-nearest neighbor hopping exceed some threshold value.) 



APPENDIX E: DERIVATION OF THE HUND'S RULE HAMILTONIAN 

The Coulomb exchange terms for the ^-states can be written in the form 28 

i a/3 aa' 

where Jh is the Hund's rule coupling. Adding TL cex to the Hamiltonian Eq. (1), the perturbation expansion in power 
of the transfer integrals t now contains a term of the order t 2 Jh/U 2 , which reads 
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(ij) Pl^(ij) crcr ' 

_ C i/3cr' C i/3o- C ] 7( T C J7cr' ~ ^C^C^ct^, Cj 7<7 / + 2cJ 7(J , C^g. C^g. Cj 7(T / . 

Taking the thermal averages using Eq. (16) we find 

(SHkk) = E E { 2A Pl {i)A 01 {j) + 8A 01 (i)A^(j) - lOA 0f3 (i)A 71 (j) 

- 2B fj . t {i) ■ B 07 (j) + AB 01 {i) ■ B 7/5 (j) - 2B p{i) ■ B 77 (j)) . 
where terms independent of the order-parameters were omitted. 



31° 



(E2) 



(E3) 



APPENDIX F: SIXTH-ORDER ANISOTROPY IN THE ORBITAL SECTOR 



At fourth-order, the terms in ai(i) and 02(1) are proportional to [a\(i) + a 2 (i)] 2 [see Eq. (41)], and there is complete 
isotropy in a\ — a 2 space. However, this isotropy must be broken in view of the special role played by the directions 
along the cubic crystal axes. This symmetry is found in the sixth-order terms, as we now show. There are several 
contributions to the free energy at sixth order in a\(i) and a 2 (i), some of which involve coupling to non-critical 
variables. To illustrate the symmetry of these terms we explicitly consider only the "direct" terms arising from Eq. 
(19), from which we have 



SF = aE TrX6 W > 

i 

where a is a numerical coefficient times kT. Thus we write 



SF = E Tr 



= a^trA 6 (i) 



(Fl) 



(F2) 



where here the trace operation, indicated by "tr," refers to a diagonal sum over the indices of the matrix A, as 
contrasted to the trace used elsewhere in this paper over the 6 i2g-states. Using Eq. (8), this yields 



01 



(*) + Via 2 (i) 



V6 



a\{i) - \Z3a 2 (i) 
V6 



+ 



-2ai(i) 



(F3) 



Now, since we are only interested in how this term affects the critical variables, we may replace VNa n (i) by a„(Q), 
which we denote a n . Then we may write 



SF 



367V 2 



10[< 



a\ + aj} 3 +a\- I5a\a 2 2 + lha\a\ - a% 



To clarify the anisotropy of this form we set a\ — r cos 6*q and a 2 = r sin 6q , in which case 

„6 



8F 



ar" 
36N 2 



[10 + cos(66> Q )] . 



(F4) 



(F5) 



This free energy has minima at the angles 9q = n/2 + nir/3, for n = 0, 1 ... 5. These correspond to a\ = —r sin(n7r/3) 
and a 2 — rcos(n7r/3). For n = 0, only a 2 is nonzero. From Eqs. (12) one sees that this corresponds to (N z (i)) = 1/3, 
and having N x (i) — N y (i) oscillate at wavevector Q with an amplitude proportional to r. By similarly analyzing 
the other minima, one concludes that these six minima correspond to the six ways one can chose indices so that 
(N a (i)) = 1/3 and (Np(i) — N 7 (i)} oscillate at wavevector Q. (There are three ways to choose a and two ways to fix 
the phase of the orbital density wave.) However, additional contributions to the free energy might make the coefficient 
of the cosine term in Eq. (F5) negative, in which case the minima occur for 8q = nir/3. Now for n = only cti 
is nonzero, and, from Eqs. (12), this corresponds to N x (i) — N y (i) = | + S(i), and N z (i) = | — 2S(i), where 8(i) 
oscillates at wavevector Q. The other minima correspond to cyclic permutations of coordinate axes consistent with 
cubic symmetry. 
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